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ABSTRACT 


Consider  a  queueing  system  in  which  customers  arrive  in 
accordance  with  a  renewal  process  having  an  unknovni  distribu¬ 
tion  F  ,  and  in  which  the  service  times  are  independent  and 
have  unknown  distribution  G  .  We  assume  that  there  are 
k(k  <_  °°)  servers. 

Let  C(t)  denote  the  number  of  customers  in  the  system 
at  time  t  .  It  is  shown  that  F  and  G  are  identifiable 
from  (C(t)  ,  t  _>  0}  if  either  G  or  F  is  continuous,  if 
F(x)  <  1  for  all  x  ,  and  if  the  number  of  busy  periods  is 
infinite  almost  surely.  Secondly,  it  is  shown  that  F  and 
G  are  identifiable  if  G  is  not  lattice  and  the  queue  size 
a.s.  converges  to  infinity.  The  second  result  is  obtained 
by  first  considering  the  related  problem:  if  N^(t)  , 
i  =  1,  ...,  k  ,  are  independent  renewal  processes  with  the 
same  distribution  G  ,  then  construct  an  a.s.  convergent 
estimate  of  G  based  only  on  {N(t)  ,  t  >_  0}  ,  where 
k 

N(t)  =  l  N  (t)  . 
i=l 


IDENTIFIAB7LITY  IN  Gl/G/k  QUEUES 
by 

Slieldon  M.  Koss 


0.  INTROIM 1CTI0N  AND  SUMMARY 

Consider  a  queueing  system  in  which  customers  arrive  in  accordance  with  a 
renewal  process  having  an  interarrival  distribution  F  ,  and  in  which  the  service 
times  of  customers  are  independent  and  have  distribution  G  .  Moreover,  suppose 
that  there  are  k(k  <_  «°)  servers  and  that  an  arriving  customer  is  immediately 
served  if  he  finds  one  of  the  servers  free,  and  if  not  then  he  joins  the  queue. 

When  a  customer  completes  service,  he  is  immediately  replaced  by  a  customer 
waiting  in  the  queue.  However,  we  do  not  assume  that  the  queue  discipline  (whether 
LIFO,  FIFO,  or  anything  else)  is  known. 

Let  C(t)  denote  the  number  of  customers  in  the  system  at  time  t  .  The 

problem  of  interest  is  to  determine  whether  F  and  G  are  identifiable  from  a 

single  sample  path  of  the  process  {C ( t )  ,  t  >_  0}  .  In  Section  1,  we.  show  that  F 

and  G  are  identifiable  if  either  G  or  F  is  continuous,  if  F(x)  <  1  for  all 

x  ,  and  if  the  number  of  busy  periods  is  a.s.  infinite.  In  Section  2,  we  show  that 

F  and  G  are  identifiable  if  G  is  not  lattice  and  if  the  queue  size  a.s. 

converges  to  infinity.  This  is  done  by  first  considering  the  following  problem 

which  is  of  independent  interest;  If  N^t)  ,  i  »  1 .  k  ,  are  independent 

renewal  processes  with  the  same  distribution  G  ,  then  construct  an  a.s.  convergent 

k 

estimate  of  G  based  only  on  the  process  (N(t)  ,  t  >_  0}  ,  where  N(t)  =  £  N,(t) 

i=l 

The  problem  of  identifying  F  and  G  has  received  some  attention  in  the 
literature.  In  [1],  Brown  considered  the  problem  of  estimating  G  from  a  single 
sample  path  of  (C(t)  ,  t  >_0)  for  the  case  M/G/°°  (that  is,  where  F  is  assumed 
exponential  and  k  *=  «)  .  By  a  clever  use  of  the  ergodic  theorem,  he  showed  that 
his  estimate  a.s.  converged  (pointwise)  uniformly  to  G  ,  which,  of  course,  shows 
that  G  is  identifiable.  Brown  also  raised  the  question  of  the  identif iabi ] ity  of 
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G  from  (C(t)  ,  t  ^  0}  for  the  GI/C/a’  system.  In  [3],  Kendall  and  Lewis  show 
that  in  a  Gl/G/“  queue  if  one  is  given  the  output  epochs  along  with  the  serial 
number  of  the  input  corresponding  to  each  output  then  the  input  distribution  is 
completely  identifiable  and  the  output  distribution  is  identifiable  up  to  a 
location  parameter. 
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1,  RESULTS  WHEN  BUSY  PERIODS  ARE  A.S.  QjjlTK 

We  say  that  a  busy  period  begins  when  a  customer  arrives  and  finds  the 
system  empty.  In  this  section  we  shall  suppose  that  a  busy  period  is  a.s.  o£ 
finite  length.  From  results  of  Kief er-Wolfowitz  [ A ,5  J  it  follows  that  this  will 
be  the  case  if  p„  <  kp„  and  F(x)  <  1  for  all  x  ,  where  p  and  pt,  are  the 
mean  service  and  interarrival  times. 

Call  a  customer  a  B-customer  if  there  are  no  other  customers  in  the  system 

when  he  arrives.  Let 

1  if  n^1  B-customer  completes  his  service  before 
a  new  customer  arrives 
0  otherwise  . 

Then,  by  the  strong  law  of  large  numbers  it  follows  that 


(1) 


I,  +  ...  +  I 
1  n 

n 


G(y)dF(y) 


Now  from  the  set  of  all  B-customcrs  let  us  consider  only  those  for  which 

there  is  an  arrival  within  a  time  x  following  their  arrival.  Call  these 

B  -customers  and  let 
x 

!1  if  B^-customer  completes  his  service  before 

a  new  customer  arrives 
0  otherwise  . 

Then  from  the  strong  law  of  large  numbers  it  follows  that  either  F(x)  =  0  and 

there  are  no  B  -customers  or  else 
x 


(2) 
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J.  +  ...  +  J 
_ n 

n 


jf  G(y)dF(y) 

n .  s .  *0 _ 

K(x)“ 


Now,  from  the  set  of  l$-cust outers  consider  those  who  finish  before  n  new 

customer  arrives,  and  let  their  service  times  be  X. ,X„ .  If  we  let  H  (x) 

n 

be  the  empirical  c.d.f.  of  X^,  ....  ,  then 


x 

/  G(y)dF(y)+C(x) (l-F(x)) 

<J)  H  (x)  a  Vs  *  H(x)  l  - - - - 

n  " 

f  G(y)dF(y) 

0 


Finally  let  ,  i  >  1  ,  be  the  successive  interarrival  times  and  let 

F  be  the  empirical  c.d.f.  of  Y. ,  ...,  Y  .  Then 
u  ‘  1  n 


(A)  Fn(x)  V*  F (x)  . 

Theorem  1: 

If 

(i)  the  number  of  busy  periods  is  a.s.  infinite 

(ii)  either  F  or  G  is  continuous 

(iii)  F (x)  <  1  for  all  x 

then  F  and  G  are  identifiable  from  {C(t)  ,  t  >  0}  . 

Proof ; 

From  the  continuity  assumption  (ii)  it  follows  that  an  arrival  and 
departure  may  not  occur  simultaneously,  and  as  a  result,  the  left  sides  of 
(1)  -  (4)  are  all  obtainable  from  the  data  (C(t)  ,  t  >  0}  .  Hence,  if 
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F(x)  +  0  then  the  right  sides  of  (1)  -  (4)  are  identifiable  and  the  result  follows 
H(x)  f  G(y)dF(y)-J  G(y)dF<y) 

since  G(x)  ■  - -0- ■  ■■■  j—— -  —  —  ■ —  .  On  the  other  hand,  if  F(x)  ■  0 

then  the  result  follows  from  (1),  (3)  and  (4). 

Q.E.D. 

Remark ; 

The  purpose  of  the  continuity  assumption  (ii)  is  that  it  enables  us  to 
determine  the  arrival  and  departure  epochs  from  the  data  { C ( t )  ,  t  >  0}  . 

If  the  experimenters  data  consists  of  the  arrival  and  departure  epochs  then  this 
condition  is  not  necessary. 


(. 


2.  THE  CASK  ku„  >  u_ 

(. _ r_ 

In  this  section  wt*  aitppose  Hull  G  is  uonlattice,  either  C  or  F  in 

continuous  and  y  *  kur  .  Pofon*  allowing  that  E  and  C!  arc  identifiable  in 
C  t  I' 

this  case  we  first  consider  the  following  related  problem. 

Related  l'roblom: 

Consider  k  independent  stationary  renewal  proe.esse3  {N^(t)],i  *  1,  k  , 

for  which  the  time  until  the  first  renewal  has  a  distribution  G  (x)  " 

e 

and  the  remaining  Jnterarrival  times  have  distribution  G  ,  where  G  is  assumed 

k 

to  be  nonlattice.  Let  N(t)  ■  £  N  (t)  .  The  problem  is  to  estimate  G  from 

i»l 

the  data  { N ( t )  ,  t  >  0}  . 

Let  E  denote  the  time  from  n  until  the  next  jump  of  the  N(l)  process, 
n 

and  let 


Lemma  .1 : 


<1 


1 

n 


=  ( 


0 


E  >  x 
n 

E  x  . 
n  - 


n-1 

l/n  l  I 
j=0  J 


Proof : 

From  known  renewal  theory  results  it  follows  that  {En  ,  n  >  0}'  is  a 

stationary  process.  Also,  by  known  renewal  results,  it  follows  that  the  distribution 

of  {E  ,  n  >  m}  given  the  values  E  ,  n  <  j  ,  converges  as  m  >  00  to  the 
n  n  - 

distribution  of  {E  ,  n  >  0}  regardless  of  the  values  E^  ,  n  <  j  .  This  implies 

that  the  tail  o-field  of  {K  ,  n  >  0}  is  independent  of  {E  ,  n  <  j}  and 

n  -  n  — 
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hence  every  tail  event  has  probability  0  or  1  which  implied  Mint  the  K 
process  is  ergodic.  Since  the  process  is  stationary  and  ergodic,  so  is  the 

I  process.  The  result  then  follows  from  the  ergodic  theorem  since 


Now,  define  T  to  be  the  time  between  the  (n  -  l)sL  and  the  n 
n 

Jump  of  the  {N(t)  ,  t  >  0}  process,  and  let 


th 


J 

n 


(1  T  >  x 
n 

0  T  <  x  . 
n  - 


Lemma  2: 


1/n  T  J,  a-V8’  (1 

i-0  1 


Proof : 


th 

Let  S  .  denote  the  time  of  the  n  event  of  the  process 

n,J 

N. (t)  ,  J  ■  1,  ...,k,n>0  .  Let  T  .  be  the  time  from  S  .  to  the  next 
J  n,j  n,j 

jump  in  the  N(t)  process.  Now  because  of  the  independence  and  stationarity 


of  the  N^(t)  processes,  1  <  i  <  k  ,  it  follows  exactly  as  in  Lemma  1  that  for 
each  j,l<j<k,{T  .  ,  n  >  1}  is  a  stationary  ergodic  process.  Letting 


1  T  >  x 

nj 

we  have  by  the  ergodic  theorem  that 

0  T  .  <  x 
nj  - 


n-1 


i  k— 1 


1/n  l  J  a-^8,  (1  -  G(x))[  /”^dx)  ,  j  =  1,2 . k  . 


i“0 


8 


The  result  row  follows  a i nee  {J^  ,  n  >  1}  Is  composed  of  the  k  subsequences 
{Jnj  .  "  2  U  ,  1  ^  J  1  k  •  Q.E.I). 


Remark : 


In  the  proof  of  l.emma  2  it  is  necessary  to  introduce  the  ^  since  the 
process  {T^  ,  n  >  1}  is  not  a  stationary  process. 

From  Lenimas  1  and  2  we  obtain 


Theorem  2: 


and  hence  C.(x)  is  identifiable  from  {N(t)  ,  t  >  0}  . 
Let  us  now  return  to  our  F/G/k  queueing  system. 


Theorem  3: 

If  G  is  nonlattice,  either  G  or  F  is  continuous,  p„  and  p  are 

G  r 

finite  and  satisfy  p  >  kp  ,  then  G  and  F  are  identifiable  from 

G  J.* 

{C(0  ,  t  >  0}  . 


Proof : 

Since  p  >  kp„  it  follows  that  with  probability  one  the  queue  size  will 
converge  to  infinity,  and  hence  after  some  finite  amount  of  time  the  output 
process  will  behave  as  the  superposition  of  k  independent  renewal  processes 
each  having  interarrival  distribution  G  .  Our  continuity  assumption  assures 
that  an  arrival  and  departure  will  not  occur  simultaneously  and  hence  departure 
epochs  will  correspond  one-to-one  with  points  of  decrease  of  C(t)  .  By  Theorem  2 


y 


it  follows  that  G  is  identifiable.  Since  points  of  increase  of  C(t) 
correspond  in  a  one-to-one  manner  with  arrival  epochs  it  follows  that  F 
identifiable  by  the  empirical  c.d.f. 

Remark; 


is 

Q .  E .D . 


As  in  the  remark  following  Theorem  1,  we  note  that  if  the  data  consists  of 
the  arrival  and  departure  epochs  then  the  continuity  assumption  may  be  dropped. 
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